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ABSTRACT 

The binary characteristics of asymptotic giant branch (AGB) stars are im- 
printed in the asymmetric patterns of their circumstellar envelopes (CSEs). We 
develop a simple method for constraining the orbital parameters of such binary 
stars from the characteristics of a spiral-like pattern observed at large distances 
from the central stars. We place constraints on the properties of AFGL 3068 (i.e., 
the masses of binary components, the viewing inclination of the orbital plane, as 
well as the orbital period, velocity, and separation). In particular, the mass of 
the companion star of AFGL 3068 is estimated to be greater than 2.6M . This 
method is applicable to other AGB stars, providing a step toward understand- 
ing the role binary stars can play in explaining the diverse patterns in observed 
CSEs. 

Subject headings: binaries: general — circumstellar matter — stars: AGB and 
post-AGB — stars: carbon — stars: late-type — stars: winds, outflows 



1. INTRODUCTION 

At the end of a star's life, low- to intermediate-mass stars return their mass to interstellar 



space by means of slow, dense winds (jHabingl Il996h during the AGB phase. The material 



expelled from these evolved stars can take the form of irregular /regular patterns and shapes 
prior to mixing with the interstellar medium. The dense CSEs often show strong asymmetries 
in the late stellar evolutionary stages, whose origin probably lies in their progenitor systems. 
Recent studies reveal that a binary system may play a crucial role in forming the various 
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shapes of CSEs (IHuggina 120071 ). However, the determination of the system properties of a 
wide binary is challenging. 

AFGL 3068 is the first AGB binary system that clearly shows a wide spiral pattern. 
Two point- like sources in nea r-infrared images have been revealed using the Keck adap- 
tive optics (IMorris et al.l 120061 ). making AFGL 3068 an evolved stellar system in which the 
binary components are indeed detected. Based on a distance of 1 kpc, the projected sepa- 
ration is 109 AU, which is relatively large to infer a direct binary interaction. Nevertheless, 
AFGL 3068 exhibits evidence of a binary interaction in view of the large scale spiral pattern 
withi n its CSE extending over 10,000 AU in dust-scattered Galactic light (IMauron fc Huggins 
20061 ). This spiral pattern provides support for t he role o f a companion star in shaping the 
CSE JSokerl Il994i iMastrodemos fc Morrisl Il999l : Eel 120071 : fedgar et al.l l2008h . n ecessitating 



an attempt to const rain binary properties from the large scale density pattern (iKiml 12011 
Kim fc Taan]l2012al Jbll 



In a parameter study of hydrod ynamical simu l ations on th e asymmetric circumstellar 
pattern due to the binary motions (IKim fc Taaml l2012bl . IKTI hereafter), we have shown 
the directional dependence of the wind flow blown out from the AGB star. It reflects the 
difference between the line-of-sight velocity of the flow and the propagation speed of the 
pattern in the plane of the sky. The pattern propagation speed is < V w > +2V p /3 throughout 
the orbital plane, but equals to < V w > in the vertical direction, resulting in an overall 
oblate shape. This shape is characterized by the aspect ratio a/b corresponding to (< V w > 
+2V p /3)/ < V w >, which appears as an ellipse with inclination. In addition, the maximum 
fluctuation of the gas velocity relative to the average wind velocity value < V w > is less than 
the orbital velocity V p of the mass-losing star. 

Based on these findings, we propose a new method to constrain the binary properties 
from four observed quantities corresponding to the projected axial ratio (a/b) pro j, the pro- 
jected binary separation (r p + r comp ) pro j, the position angle of the star from the major axis of 
the elliptical shape (PA), and the arm spacing Ar arm . Here, r p and r comp correspond to the 
separation of the AGB star and its companion from the center of mass of the system. For 
the first application of this method, we place constraints on AFGL 3068, as an extreme car- 



bon star (Lebofskv &: Riekel 1977 ; Jones et al. 1978 ; Silverglate et al. 1979 ; Jewell fc Snyder 



19821 : IVolk et al.lll992l ; Izhang et al.ll2009l ). at the tip of the AGB. Figured] shows the HST 



image of AFGL 3068, overlaid with our binary model with the mass ratio of 1.13. 
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Fig. 1. — (a) Image of AFGL 3068 in the wide-V F606W filter of the Hubble Space Telescope- 
ACS camera with the field size of 25% x 25% (Credit: Mauron & Huggins, A&A, 452, 257, 
2006, reproduced with permission © ESO.). North is up, and east is left. The bright 
straight line refers to a diffraction spike of a st ar in the field of view of the original image, 
(b) Same image with an Archimedes spiral fit ( iMauron fc Hugginsl 120061 1 . overlaid with our 
three-dimensional hydrodynamic model. The observed pattern is more elongated than an 
Archimedes spiral having a constant spacing, particularly in the northwest direction. The 
elongated shape of the spiral pattern is described by the ellipses with the axial ratio of 1.1 
corresponding to the eccentricity of 0.417 (dashed lines). The ellipses have the major axis 
(the intersection of the plane of the sky; line of nodes) at 70° countercloc kwise from east- 
west, where the binary stars are aligned at the center (Morris et al.l 120061 1 . and the minor 
axis is perpendicular to the major axis (black thin lines). The difference in the semi- major 
radius between the dashed lines is 2350 AU at a distance of lkpc. (c) Density distribution 
at the mid-plane in the line of sight of a hydrodynamic model for binary stars composed of 
a mass-losing star (3.5M ) and a main-sequence star (3.1M Q ) with the orbital separation 
of 166 AU. The corresponding orbital velocities are 2.8 km s -1 and 3.2kms _1 , respectively. 
The wind velocity varies about a mean value of 11.5 km s -1 . An adiabatic index 7 of 1.4 
is employed to produce a te mperature profile o f the form r -0 83 as obtained from the CO 
molecule excitation analysis ([Woods et al.ll2003l ). The orbital plane is inclined by 50° from 
the plane of the sky. (d) Column density distribution. 
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Fig. 2. — Schematic diagram for the process of deriving binary properties (open boxes) from 
four observables (pink boxes), based on a parameter study using hydrodynamic simulations 
(IKTl ). The speeds of orbital motion V p and wind motion < V w > of the mass- losing star 
are considered as the model parameters. The four required observables correspond to the 
projected major-minor axial ratio of the ellipse characterizing the overall shape of the pattern 
(a/b) proj, the projected orbital separation (r p + r com p) pro: ,-, the position angle of the stars with 
respect to the major axis (line of nodes), and the spiral arm spacing Ar arm . The result of 
the analysis leads to the masses M p and M comp , orbital radii r p and r comp for the mass- losing 
star and the companion, respectively, the orbital period T p , and the inclination angle of the 
orbital plane i. 
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2. NEW METHOD 

Figured is a schematic diagram summarizing the procedure of constraining the binary 
stellar and orbital properties, which takes V p and < V w > as two free parameters. Consider 
an oblate spheroid given by the relation (x/a) 2 + (y/a) 2 + (z/b) 2 = 1 to describe the overall 
shape of the circumstellar pattern in a three-dimensional space. The determination of the 
projection is key in connecting the observed features to the binary properties. The projected 
shape of the oblate spheroid is obtained by coordinate transformation to the coordinate 
rotated with the inclination angle i. Defining the the line of nodes ( 

and Uproj = ycosz+zsinz) yields the elliptical shape in the projected plane as (x pro j/a pro j) 2 + 
(Vproj/bproj) 2 = 1 where a pro j = a and b~ 2 Q - = (cosi/a) 2 + (sini/6) 2 , which simplifies to 

2 (a/&fro 3 -i m 
sm 1 ~ ( a/6 )2 _ x • w 

Hence, the orbital inclination can be estimated from the observed (projected) axial ratio, 
given the relation between the aspect ratio of the oblate spheroid and the orb ital-t o-wind 



velocity ratio, a/b = 1 + 2V p /(3 < V w >), as found from our parameter study in lKT 



Similarly, simple geometry yields the orbital separation of the binary stars (r p + r comp ) 
as a function of the inclination i. Figured illustrates the inclined orbit with respect to the 
plane of the sky, and we obtain the projection effect for the orbital radius of a star to be 
given by r PtPro j = r p (l — sin 2 cf) p sin 2 i) 1 ^ 2 and the same for the companion's orbital radius 
r comp . Therefore, the true orbital separation can be obtained using the observed separation 
(r p + r comp ) pro j and position angle PA of the stars from 

[Tp + ?"comp )proj /„\ 
Vp + r comp — JZ . 2 7 • 2 -M /9 V / 

(1 — sm <p p sm i) v l l 
where tan0 p = tan(PA)/ cosi, 

and the inclination derived from Equation (JT]). 

On the other hand, the arm spacing Ar arm of the spiral pattern in a binary system is 
determined by the product of the binary orbital period and the pattern propagation speed 
in the orbital plane: 

/ 2 \ 2irr 

Ar arm = < V w > +-V P x — (3) 



3 7 v P 



This equation shows that the orbital radius r p can be obtained from the pattern spacing 
Ar arm along the observed longest axis, which refers to a section of the orbital plane, for a 
given velocity ratio V p j < V w >. Using r p + r comp and r p from Equations (j5J) and (jSJ), the 
binary mass ratio q = M p /M comp corresponding to r comp /r p is determined. 



Fig. 3. — Sketch illustrating the binary orbital plane inclined from the plane of the sky with 
the inclination i. The separation of the mass- losing star and companion from the common 
center of mass are r p and r comp , respectively. The longitudinal angle from the line of nodes 
is (j) p in the orbital plane, and the projected angle in the plane of the sky is defined by PA. 
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Finally, we employ the formula for the force balance in a binary system, 

T/2 CM 

v p UJ ' 1 comp / j \ 

'"p V p "I - '"comp) 

and replace the quantities related to the companion by the quantities of the mass-losing star 
using the definition of the q factor, resulting in M p = q(l + q) 2 V p r p /G. The companion mass 
is M p /q. 



3. CASE STUDY OF AFGL 3068 

3.1. Observed Properties 

We first find the best fit of an elliptical shape with the dust scattered light pattern of 
AFGL 3068 (Figure[T^b)). The modeled ellipse is characterized by an axial ratio (a/b) pro j of 
1.1 with the major axis in 70° from the west (i.e., right). Identifying the major axis as the 
line of nodes of the inclined orbital plane for the reference direction, the position angle (PA) 
of two point-like sources arranged to the east- we st direction is 70° . The observed separation 



between two sources (r p + r comp ) pro j is reported (IMorris et al.l 120061 ) to be 109 AU (<i/lkpc), 
and the interval between ellipses along the major axis, Ar arm , is 2350 AU (d/1 kpc), where d 
refers to the object distance. The slight deviation from a very good match of the elliptical 
shape with most parts of the spiral pattern implies that partial circular (or elliptical) rings 
observed in other AGB CSEs could also contribute to the spiral pattern driven by a binary 
motion. 



3.2. Parameter Space 



For AFGL 3068, we construct a parameter space, in which V p and < V w > are spanned, 
and follow the procedure summarized in Figured to calculate the corresponding orbital 
properties of the binary system. FigureH] shows the result, restricting the parameter space 
that satisfies the conditions for AFGL 3068 as a carbon-rich AGB star in a binary system. 
First, the gas velocity cannot exceed < V w > +V P in any direction (IKTI ). which yields 



< V w > +V P > 14 km s 1 with the line- of-sight velocity of 14kms 1 defined from the half 



width of spectral lines (ILoup et al.lll993l ). Second, the mass of a carbon star such as AFGL 
3068 is theoretically constrained to be less than 4M Q since more massive stars (4 to 8M ) 
are sufficiently hot at the base of their convective envelopes to activate t he so-called "hot 



bottom burning" process, which significantly burns the envelope carbon (IBoothroyd et al. 




V p (km/s) 



Fig. 4. — Model parameter space for the speeds of orbital motion V p and outflowing wind < 
V w > of the ma ss-losing star in A FGL 3068. An expansion velocity obtained from molecular 
line emissions ( jLoup et al.lll993l ) provides a constraint, < V w > +V P > 14kms~ 1 (dotted 
line), and the mass of the AGB carbon star (IBoothroyd et al.lll993l . 119951 ) can be 0.6-4M© 
(black solid line), confining the parameters for AFGL 3068 into the gray-colored region. 
Considering the carbon star to be more massive than its companion, the parameter space 
for AFGL 3068 is further restricted to above the red line for the mass ratio of 1. The + 
sign corresponds to the parameters for which a hydrodynamic simulation was carried out 
for comparison to observations (Figured]). Top left. The mass ratio between the carbon 
star and companion of 0, 1, and 2 (red solid lines) and 0.2, 0.4, 0.6, and 0.8 from bottom 
(red dashed lines). Top right. The mass of the companion with an 0.5M Q -interval (yellow 
dashed lines). Bottom left. The inclination angle of the orbital plane with respect to the 
plane of the sky with a ten-degree interval (green dashed lines). Bottom right. The orbital 
period of binary stars with a 50- year interval (blue dashed lines). 



-9 - 



19931 . Il995[ ). Third, w e take the lower lim it of the carbon star as O.6M , which is the mean 
mass of white dwarfs (IMadej et al.l I2004J ) . 



These three constraints place the AFGL 3068 system within the gray-colored area in 
FigureHl As a result, the mass of the companion is constrained to be greater than 2.6M 
with the binary mass ratio q = M p /M comp greater than 0.2. In addition, the inclination angle 
of the system is constrained to lie within approximately 40°-50°. Finally, the orbital period 
is found to be less than 870 years, and the orbital speed of the carbon star in the range of 
2.5-4 km s -1 . Assuming an accuracy of the distance measurement to be 20%, the companion 
mass is est imated to be greater than 2.1 to 3.2M . If the line-of-sight velocity of the wind is 
16kms _1 ( IDe Beck et al.ll2010f ). the estimated companion mass is larger than 2.8 or 4.2M 
at distance of 0.8 or 1.2 kpc. The combination of a high wind velocity (16kms _1 ) and large 
distance (1.2 kpc) results in a large companion mass, leading to an implausible mass >4M , 
for which the companion should have evolved more rapidly than the carbon star. 

As more massive stars evolve more rapidly, the carbon star is likely to be more massive 
than its companion (q > 1), although the mass ratio q could be less than unity because the 



carbon star is in t he phase of losing mass at a high rate (> 10 5 M^yr 1 ) ( Jones et al. 19781: 



2008; 



| Ladjal et al. 



Volk et all 119921 ; iWoods et al.l 120031 : lLoup et al.l 119931 : iRamstedt et al 
2010). In the case that the total mass lost from the carbon star is not too extensive, the 



parameters can be further constrained into a small parameter space, indicating individual 
stellar masses of 3-4M , speed of 2.5-3.5 km s -1 with the orbital period of 700-850 years in 
the plane inclined by approximately 50°. 



3.3. Hydrodynamic Model 

Within the limited area in the parameter space where AFGL 3068 can lie, we select a 
point (plus sign in FigureH]) to compare the corresponding hydrodynamic model with the 
observed pattern in the CSE of AFGL 3068. FigureH] indeed shows that the model shows 
good agreement with the scattered light image of AF GL 3068. In particular, t he Archimedes 



spiral characterized by having a constant spacing ( iMauron fc Hugginsl |2006[ ) passes inside 
the intensity peak positions in the northerly direction, while the spiral of our model better 
traces the positions of the pattern. The blurred patterns seen in the northwest and south 
closely resemble the knotty structures appearing along the line of nodes, which form from 



the overlap of the two structures generated by the orbital motions of individual stars (JKTJ). 



The modeled binary system is characterized by a mass-losing star of mass M p = 3.5M , 
velocity V p = 2.8 km s -1 , and orbital radius r p = 78 AU, and a companion of mass M comp = 
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3.2 km s 1 , and orbital radius 



comp 



88 AU. The wind velocity 



3.1M Q , velocity V comp 

V w was set to be 10.9 km s -1 at 30 AU distance from the mass-losing star. The gas flow 
evolves to a quasi-steady state fluctuating about 11.5 km s -1 with a standard deviation of 
0.67 km s -1 . To achieve a te mperature profile o f r -0 83 as determined from an analysis of 
the CO molecule excitation (IWoods et al.l 120031 ) for a density profile of r~ 2 , the adiabatic 
index 7 of the wind material is set to be 1.4 with the initial temperature of 1540 K (/x/2) at 
30 AU distance from the mass-losing star, where /1 indicates the molecular weight. Figure[5] 
shows that the resulting temperature profile indeed follows the expected power law with the 
index of —0.83 and has an average value of approximately 60K(/i/2). Additional heating 
in the post-shock region and cooling via molecular li ne transitions are n ot considered as the 
timescales in the CSEs of AGB stars are comparable (jEdgar et al.ll2008l ). The density power 
law index is —2.2 in the resulting configuration. The volume and column density contrasts 
correspond to a factor of 10-30 and 2-5, respectively. We note that the density contrast 
in the model with 7 = 1.4 is approximately constant, while th e con trast in an isothermal 
model decreases rapidly with distance from the system center ( iKTi ). The dependencies of 
the overall shape and morphology of the spiral pattern on the orbital-to-wind yelocity ratio 
are in sensitive to the value of 7. The mass loss rate taken as 4 x 10 -5 M & yr _1 (jOlivier et al. 



200ll ) provides the scale for the magnitude of the density, but not the density contrast. The 



resulting column density toward the star is of order of 10 cm , which varies by a factor of 
a few with the inclination angle of the system. 

The observed pattern thickness of AFGL 3068 is thinner than the arm thickness of the 
modeled column density distribution, and is rather comparable to the density distribution 
in the mid-plane of the sky. Because of the curved geometry of the spiral-shell pattern, the 
shocked edges of the pattern at different layers are projected into different positions in the 
plane of the sky, resulting in a broader arm for the column density. It may imply that the 
observed dust scattered light is optically thick so that the observed image refers to only part 
of the pattern in the line of sight. This is co nsistent with the h igh optical depth derived 
from the spectral energy distribution fitting ( ILadjal et al.l 120101 ) . The observed pattern 
thickness may also depend on the observed waveband. In particular, the size distribution 
of the dust grains may be affected, resulting from coagulation or shattering via grain-grain 
collisions immedia tely after encountering the shock with a velocity dispersion of a few kms -1 
( llones et al.lll996l ). As a consequence, the variation of the size distribution may affect the 
intensity across the shock pattern at the observed wavelength. Small grains with a size less 
than 0.045 /zm tend to strongly couple with the gas, while larger gra ins are characterized b y 
a different distribution as they are more resistant to the drag force (Ivan Marie et al.l 1201 ll ). 
The probable large population of small grains near the shock increases the scattering cross 
section with the large surface-to- volume ratio of the small grains, resulting in an increase of 
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Fig. 5. — Temperature profile in the outflowing CSE surrounding the binary system 
along two arbitrary-chosen directions in the orbital plane (thick and thin solid lines) in 
the hydrodynamical model. The peak temperature within the spiral arm has the value 
2.8 x 10 5 K (/z/2)(r/AU)~ ' 83 (thick dashed line), while the inter-arm temperature approxi- 
mates to 1.5 x 10 4 K (/i/2)(r/AU)~ a83 (thin dashed line). The temperature along the orbital 
axis (dot-dashed line) and the value averaged over three dimensional space (dotted line) are 
consistent with the result of a radiative transfer modeling for the observed CO line emission 
JWoods et al.ll2003h . 
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the scattered light intensity at the corresponding wavelengths. 



Our model assumes circular motion of the binary stars. However, previous calculations 
for eccentric orbits do no t show a noticeable change in the overall shape outside the second 
turn of the spiral pattern (jHell2007l ; iRaga et al.ll201ll ). suggesting that our model for the large 
scale shape of the pattern driven by circular orbits is a good approximation even for the case 
of non-circular orbits. Although a highly eccentric orbit (e > 0.8) causing substructure in 
the inter-arm region (prominent between the innermost turns) may be responsible for the 
extra structure found in the observed image between first and second turns in the northwest 
and southeast, a circular orbit model also exhibits a weak density enhancement close to the 
observed substructure. This substructure can also possibly be affected by local phenomena 
around the companion (e.g., the formation of an accretion disk and radiation pressure effects). 



4. Conclusion 

We have developed a new method to estimate the stellar and orbital characteristics of 
binary systems with an AGB component using information gleaned from the properties of 
a spiral pattern in the CSE. As a first application, we have provided constraints on the 
inclination angle, orbital period, and companion mass in the binary system AFGL 3068 
based on the pattern shape of AFGL 3068 seen in dust scattered light. With molecular line 
observations, additional kinematic information on the velocity structure of the pattern will 
provide further constraints on the system properties. 

The potential for estimating system parameters with the method proposed in this paper 
is promising and can be useful in providing fundamental input about systems where asym- 
metries in CSEs are observed. As such one may be able to use these observational data to 
place constraints on the input parameters and prescriptions used for the evolution of syst ems 



in population synthesis studies of detached binary systems (see iPolitano &: Taamll201lh . It 
is likely that the application of this method to other systems, which we can anticipate will 
be found with very sensitive and high resolution observations using e.g., the Atacama Large 
Millimeter /submillimeter Array, will provide a greater statistical basis for determining the 
incidence of such evolved binary systems in the total stellar binary population. 
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